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We revisit the cosmological evolution of domain wall networks, taking advantage of recent im- 
provements in computing power. We carry out high-resolution field theory simulations in two, three 
and four spatial dimensions to study the effects of dimensionality and damping on the evolution 
of the network. Our results are consistent with the expected scale-invariant evolution of the net- 
work, which suggests that previous hints of deviations from this behavior may have been due to 
the limited dynamical range of those simulations. We also use the results of very large (1024^) sim- 
ulations in three cosmological epochs to provide a calibration for the velocity-dependent one-scale 
model for domain walls: we numerically determine the two free model parameters to have the values 
c„ = 0.5 ± 0.2 and fc» = 1.1 ± 0.3. 



I. INTRODUCTION 

Phase transitions that are thought to have happened 
in the early universe have a number of inevitable conse- 
quences, the most interesting of which is the formation 
of topological defects [ll, [1| . The literature on the sub- 
ject has (for good reasons) focused on cosmic strings, 
but other defects can be of interest too. Domain walls, 
being the simplest defect (they can be described by a 
single scalar field) provide a simple testbed where one 
can study how several physical mechanisms influence de- 
fect evolution, and this knowledge can then be applied to 
other defects. This is despite the fact that the observa- 
tional roles of domain walls are very tightly constrained: 
the Zel'dovich bound Q rules them out if their symme- 
try breaking scale is > IMeV, and the bound is even 
tighter for wall networks with junctions 

Here we take advantage of the continuous improve- 
ments in computing power to carry out a large set of 
high-resolution simulations of domain walls in two, three 
and four spatial dimensions, usiiig the standard Press- 
Rydcn-Spergel (PRS) algorithm [5:1]. While the 3D simu- 
lations are of obvious cosmological relevance, the 2D ones 
have the advantage of allowing for a much larger dynam- 
ical range, and the 4D ones may be relevant to some 
brane world scenarios [^-Q. We'll only consider the sim- 
plest domain wall model, thus neglecting scenarios where 
domain walls have junctions. 

One can also describe the broad macroscopic proper- 
ties of these networks by an analytic model, in the same 
spirit of the model of Martins and Shellard for cosmic 
strings The large-scale features of the network 

are described by a lengthscale (or correlation length) L 
and a microscopically averaged (root-mean-squared) ve- 
locity V. Although this has the advantages of tractabil- 
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ity and conceptual simplicity, these come with a price: 
in going from the microphysics to the averaged evolution 
equations one is forced to introduce phenomenological 
parameters which parametrize our ignorance about cer- 
tain dynamical mechanisms, and the only way to deter- 
mine the correct values of these parameters is by referring 
to numerical simulations to calibrate them. 

The evolution of domain walls networks has been pre- 
viously studied, numerically, by a number of different 
authors [1, fl3l - lT9l |. who usually found some hints for 
deviations to the scale-invariant evolution. Since there 
are good reasons to expect this scale-invariant attractor 
18j, i20|, one may ask whether these deviations point to 
the presence of physical mechanisms not accounted for in 
analytic descriptions or if they are simply a consequence 
of the limited dynamical range of numerical simulations. 
While the present work cannot completely resolve this 
issue, we believe that it provides support for the second 
alternative. 

In what follows we briefly describe the PRS algorithm 
used in numerical simulations, as well as the analytic 
model we will use to describe the domain wall networks. 
We then proceed to present the main results of our simu- 
lations, and finally bring the two approaches together by 
comparing the simulation results to the analytic model 
predictions, thereby providing a calibration of the free 
parameters of the latter. We'll conclude with brief com- 
ments on the cosmological implications of our results and 
on possible follow-up work. Throughout the paper we 
shall use fundamental units, in which c = h = 1. 



II. DOMAIN WALLS AND THE PRS 
ALGORITHM 

We'll be interested in domain wall networks in flat ho- 
mogeneous and isotropic Friedmann-Robertson- Walker 
(FRW) universes. A scalar field (j) with Lagrangian den- 
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sity 



(1) 



where we wih take V{4)) to be a 0^ potential with two 
degenerate minima, such as 



(2) 



will have domain wall solutions. By the standard vari- 
ational methods we obtain the field equation of motion 
(written in terms of physical time t) 



dc^ 2 dv 



(3) 



where V is the Laplacian in physical coordinates, H — 
a~^{da/dt) is the Hubble parameter and a is the scale fac- 
tor, which we assume to vary as a oc t^. In what follows 
we will be interested in comparing the network evolution 
in several different cosmological epochs, in particular the 
radiation era (A = 1/2), the matter era (A = 2/3) and a 
fast expansion era (A = 4/5). 

We then apply the procedure of Press, Ryden and 
Spergel 's'l, modifying the equations of motion in such 
a way that the thickness of the domain walls is fixed in 
co-moving coordinates. One expects that this will have 
a small impact on the large scale dynamics of the do- 
main walls, since a wall's integrated surface density (and 
surface tension) are independent of its thickness. In par- 
ticular, this assumption should not affect the presence 
or absence of a scaling solution @ , provided one uses a 
minimum thickness jl7| — we will briefly revisit this issue 
below. (For a detailed discussion of analogous issues in 
the context of cosmic strings see . ) 

In the PRS method, equation ^ becomes: 



drj'^ \d\nri J drj 



(4) 



where rj is the conformal time and a and /3 are constants: 
/3 = is used in order to have constant co-moving thick- 
ness and a = 3 is chosen (in 3D, see [l^j for an argument 
in other dimensions) to require that the momentum con- 
servation law of the wall evolution in an expanding uni- 
verse is maintained ^3]. In fact we have simulated net- 
works with various values of the damping coefficient a. 

Equation (g]) is then integrated using a standard finite- 
difference scheme. We assume the initial value of (p to be 
a random variable between — 0o a-nd +(I)q and the initial 
value of d(f)/drj to be zero. This will lead to large en- 
ergy gradients in the early timesteps of the simulation, 
and therefore the network will need some time (which is 
proportional to the wall thickness) to wash away these 
initial conditions. The conformal time evolution of the 
co-moving correlation length of the network (specif- 
ically A/V oc ^J^, A being the comoving area of the 
walls) and the wall velocities (specifically 7W, where 7 is 



the Lorentz factor) are directly measured from the sim- 
ulations, using techniques previously described in [l8l |. 
However we now used a newly parallelized version of the 
code, optimized for the Altix UVIOOO architecture of the 
COSMOS Consortium's Universe supercomputer. 



III. THE ONE-SCALE MODEL 

The way to analytically model defect networks is to 
start from their microscopic equations of motion (the 
Nambu-Goto equations, in the case of strings) and carry 
out a statistical average, under the assumption that the 
defects are randomly distributed at large enough scales. 
This leads to a macroscopic energy evolution equation 
(which one can equivalently write as an equation for the 
network's correlation length or for another suitable char- 
acteristic lengthscale) and an equation for the network's 
root-mean squared (RMS) velocity. 

These evolution equations provide a 'thermodynami- 
cal' description of the network, in the same sense that 
the microscopic equations provide a statistical physics 
one. The more subtle part of this procedure is that of 
the addition of terms in these equations to account for 
defect interactions and energy losses. Such terms must 
be added in a phenomenological way, and for their cali- 
bration one must resort to numerical simulations. 

For cosmic strings, this procedure leads to the so- 
called velocity-dependent one-scale (VOS) model [lol - [l^ . 
which has been well-tested against simulations and is 
used for predicting CMB signals of string networks. One 
can follow an analogous procedure both for the case of 
monopoles and for domain walls. This latter case 
has been described in |_18], which also provided a sim- 
ple (qualitative) calibration; later in this article we will 
revisit this issue and provide a more quantitative one. 

The evolution equation for the characteristic wall 
lengthscale L (which is related to the wall density pw 
via L = cr/pu,, where tr is the domain wall energy per 
unit area) and their RMS velocity w, are as follows 



dL 
Hit 



(1 + ?,v'^)HL + c^v 



dv 
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(5) 
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the latter equation was first obtained, using a different 
approximation, in [23 |. Here and A;^ are the free pa- 
rameters, to be calibrated using simulations: the former 
quantifies energy losses, while the latter quantifies the 
(curvature-related) forces acting on the walls. At least to 
a first approximation, these are expected to be constant. 
Moreover, in the context of the VOS model the charac- 
teristic lengthscale L can further be identified with the 
physical correlation length £,phys- The comoving version 
of this was defined in the previous section; the two are 
related through 



(7) 
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and we are therefore assuming that (,phys = L. We wih 
use the two terms interchangeably for the rest of this 
paper. Note that the scaling exponents of the two cor- 
relation lengths, relative to their respective times, are 
different: if 



then 



^phys CX ^ ^ ^7 



(8) 
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for an expansion rate A defined as before. 

If one neglects the effect of the energy density in the do- 
main walls on the background (specifically, on the Fried- 
mann equations) — which is the relevant context for our 
numerical simulations — it is easy to see that, just as for 
cosmic string networks, the attractor solution to the evo- 
lution equations (I5I6P corresponds to a linear scaling so- 
lution 

L = et , V = const . (10) 

Assuming that the scale factor behaves as a (x the 
detailed form of the above linear scaling constants is 

2 ^wij^w 



3A(1 - A) 

1 — A fcu) 
3 A /l7m -\- c 



(11) 



(12) 



As in the case of cosmic strings [25| , an energy loss mech- 
anism (that is, a non-zero c^) may not be needed in or- 
der to have linear scaling. In the c^i limit one finds 
that for a > 1/4 a linear scaling solution is always possi- 
ble. Therefore, a linear scaling solution for domain walls 
can always exist in both the matter and the radiation 
eras, which shows that having non-standard (that is non- 
intercommuting) domain walls is by no means sufficient 
to ensure a frustrated wall network. 

In passing, we note that the cosmological linear scal- 
ing solutions for walls imply that the wall density grows 
relative to the background density, and will eventually 
become dominant (unless some mechanism like a subse- 
quent phase transition were to make it decay and dis- 
appear). This is ultimately the reason for the Zel'dovich 
bound W\- Moreover, since the wall density grows relative 
to the background, it must be included in the Einstein 
equations. In this case (further discussed in [3, [13) 
domain wall network will become frozen in comoving co- 
ordinates with L (X a and the scale factor growing as 
a (X t^. Notice that this solution does not depend on 
Cw — domain wall interactions play no role here since the 
walls are effectively frozen. 



IV. BOX SIZE, DAMPING AND WALL 
THICKNESS EFFECTS 



We have started by running several sets of simulations 
in two, three and four spatial dimensions, varying the box 




FIG. 1: The evolution of the domain wall density (p^ ~ A/V) 
and velocity {(■yv)'^) as a function of conformal time, for the 
3D 512'^ boxes of Table U corresponding to different amounts 
of damping, a = 2,3, 4. The plotted curves are the aver- 
age of 25 different simulations with different (random) initial 
conditions. 



size, the damping parameter a and the wall thickness 
(the typical number of points describing each wall, de- 
noted Wo), as further detailed below. These — and indeed 
all the simulations described in this paper — were started 
at a conformal time 770 = 1 and evolved in timesteps 
A77 = 0.25770 until a conformal time equal to half the 
box size, and each set of simulations contains 25 runs 
with different (random) initial conditions. Unless other- 
wise stated the quoted errors are statistical errors in the 
ensemble of 25 runs. 

Our main concern here is with a diagnostic for scaling. 
We looked for the best fit to the power laws 
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Box Size 


Q 


Fit range (rj) 


M 


V 




71) 


4096^ 


2.0 


41-2048 


-0.95 ± 0.04 


-0.00003 ± 0.00002 


0.61 


0.37 


4096^ 


3.0 


23.5-2048 


-0.97 ±0.05 


-0.00002 ± 0.00002 


0.58 


0.34 


4096^ 


4.0 


21-2048 


-0.96 ±0.04 


-0.00001 ± 0.00002 


0.55 


0.33 


8192^ 


2.0 


12.25-4096 


-0.98 ± 0.03 


-0.000011 ± 0.000008 


0.91 


0.40 


8192^ 


3.0 


21-4096 


-0.97 ±0.05 


-0.000011 ± 0.000006 


0.82 


0.32 


8192^ 


4.0 


21-4096 


-0.96 ± 0.04 


-0.000007 ± 0.000005 


0.76 


0.29 


256^ 


2.0 


38.5-128 


-0.94 ±0.09 


-0.0000 ± 0.0005 


0.56 


0.44 


256^ 


3.0 


26-128 


-0.92 ± 0.05 


-0.0000 ± 0.0003 


0.52 


0.37 


256^ 


4.0 


26-128 


-0.94 ±0.04 


-0.0000 ± 0.0002 


0.49 


0.33 


512^ 


2.0 


38.5-256 


-0.96 ± 0.05 


-0.0001 ±0.0001 


0.57 


0.43 


512^ 


3.0 


26-256 


-0.97 ±0.04 


-0.0000 ± 0.0001 


0.54 


0.37 


512^ 


4.0 


26-256 


-0.96 ± 0.03 


0.0000 ± 0.0001 


0.50 


0.34 


64* 


2.0 


21-32 


-1.0 ±0.1 


0.001 ± 0.002 


0.53 


0.47 


64* 


3.0 


21-32 


-1.0 ±0.1 


0.000 ± 0.002 


0.49 


0.35 


64* 


4.0 


21-32 


-1.0 ±0.1 


0.001 ± 0.001 


0.49 


0.35 


128* 


2.0 


38.5-64 


-1.0 ±0.2 


0.0000 ± 0.0009 


0.60 


0.47 


128* 


3.0 


21-64 


-0.9 ±0.1 


0.0002 ± 0.0005 


0.54 


0.42 


128* 


4.0 


23.5-64 


-0.9 ±0.1 


0.0003 ± 0.0001 


0.54 


0.40 



TABLE I: The measured scaling exponents \i and v (with one-sigma statistical uncertainties) for 2_D, 3D and 4_D numerical 
simulations of domain wall networks with different box sizes and damping terms (a). The third column shows the range of the 
part of the simulation that was actually used in order to fit for the scaling exponent. The last two columns show the directly 
measured asymptotic values of (,clr\ and 711, which can be related to the macroscopic parameters of the analytic model. All 
simulations were done in a matter-dominated era (A = 2/3),with /3 = and a constant wall thickness Wo = 10. Each value of 
the scaling exponent is obtained by averaging over 25 simulations with different random initial conditions. 



for a scale-invariant behavior, we should have /x = — 1 
and V = The behavior of the scaling exponent for the 
network's kinetic energy (or velocity), has not been 
explored in detail in previous work. For the scaling of 
the correlation length, it has been suggested that this 
can be slightly larger than /i = — 1, which would imply 
that the network is not evolving as fast as is allowed by 
causality. This would presumably suggest the presence of 
physical mechanisms influencing the network that were 
not being accounted for. 

One must be careful to fit only the reliable dynamical 
range of each simulation. As previously pointed out, the 
dynamics at the beginning of the simulation will be domi- 
nated by our choice of initial conditions will influence the 
evolution. Moreover, for the walls to be sufficiently well 
defined (which is certainly helpful in accurately measur- 
ing walls areas and velocities) the co-moving correlation 
length should be significantly larger than the wall thick- 
ness. Since we end all the simulations when the horizon 
becomes half the box size, the periodic boundary condi- 
tions should have no influence on our results. Our choice 
of the reliable period for the fit was done by inspection 
of each set of simulations, using these criteria 

Table |T] shows the compared results of matter-era sim- 
ulations where we varied the box size and the amount of 
damping (parametrized by a) while keeping a constant 
wall thickness as in the PRS prescription. The key re- 



sult here is that all the scaling exponents are consistent 
with the presence of a scale-invariant evolution of the 
network, corresponding to /i = — 1 and = 0. The con- 
stancy of the wall velocities (which had not been studied 
in quantitative detail by previous authors) is particularly 
noteworthy. 

We see no hints of the deviations from this scaling 
behavior that have been discussed by previous authors, 
particularly Q and [13] • (On the other hand, [l^ finds 
a possible deviation in the largest of their simulations, 
while the smaller ones are consistent with scaling.) We 
therefore suggest that such hints may have been due to 
the limited dynamical range of simulations. Interestingly, 
one can observe (particularly in the 3D case) that dou- 
bling the box size tends to bring the scaling exponent /i 
closer to —1, although it has to be said that this may in 
part be due to the (somewhat arbitrary) choice of the dy- 
namical range that is used for the fit. We will comment 
of this point again in the next section. 

In addition to showing the calculated scaling expo- 
nents /I and the table also shows the directly measured 
asymptotic values of ^c/?y and 71;, which can be related to 
the macroscopic parameters of the analytic model. These 
are calculated from the last few timesteps of each simu- 
lation, on the assumption that by then the network has 
reached scaling — note that our measured scaling expo- 
nents are consistent with this assumption. Given the sim- 
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FIG. 2: Same as Fig. [2 for the 2D 8192^ simulations. 



plicity of this method we have chosen not to present the 
(statistical) error bars on these numbers, which should 
therefore be seen as qualitative indicators of the proper- 
ties of the network. (The errors we find are typically at 
the ten to twenty percent level, with some dependence 
on the box size.) A more detailed analysis of the scaling 
properties will be discussed in the next section. 

These results also confirm the expectation that increas- 
ing the amount of damping leads to slower walls. One of 
the consequences of this is that reduces the rate of wall in- 
tercommutings and formation of closed walls. Since this 
is a key energy loss mechanism form the wall network, 
it follows that a larger damping also leads to a higher 
wall density, corresponding to a smaller wall separation 
or correlation length. 

Figs. [T] and [2] show two examples for 2D and 3D sim- 
ulations. The fact that the 2D plots show much larger 
fluctuations is due to the fact that the precision of our 
algorithm for measuring domain wall areas and velocities 
is somewhat dependent on dimensionality and box size. 
We have chosen to always use the same algorithm, de- 




_g I 1 1 1 1 1 1 

D 1 2 3 4 5 6 

Ir.(Tl) 




-8 I 1 1 1 1 1 1 

D 1 2 3 4 5 6 

Ir.(Tl) 

FIG. 3: The evolution of the domain wall density (p„ — A/V) 
and velocity {{"/v)'^) as a function of conformal time, for the 
radiation era (A = 1/2) boxes of Table UT] corresponding 
to different values of the (constant) wall thickness, Wo = 
10, 50, 100. The plotted curves are the average of 25 different 
simulations with different (random) initial conditions. 

scribed in [l8l |. rather than individually optimizing it for 
different boxes. 

Table inland Figs. |3] and |4] show the compared results 
of radiation-era (A — 1/2) and fast expansion (A = 4/5) 
simulations, applying the PRS procedure to 512"^ boxes 
and varying the thickness of the domain walls. As ex- 
pected the relaxation time of the network is directly 
proportional to the wall thickness, being approximately 
given by the light crossing time of the walls. 

Mindful of the fact that in the case of the thicker walls 
the network is just reaching the scaling solution at the 
end of the simulation, the late-time values of wall densi- 
ties and velocities are consistent with each other within 
the numerical uncertainties. Note, in particular, that 
the velocities are remarkably similar; the differences in 
the wall densities may, at least in part, be due to the 
algorithm being used to identify the walls. 
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A 


Wo 


Fit range [rf) 




V 




711 


1/2 


10 


38.5-256 


-0.99 ± 0.05 


-0.0001 ± 0.0002 


0.60 


0.46 


1/2 


50 


188.5-256 


-1.0 ±0.5 


0.0001 ± 0.0008 


1.00 


0.53 


1/2 


100 


188.5-256 


-1.0 ±0.4 


0.0003 ± 0.0007 


0.90 


0.52 


4/5 


10 


21-256 


-0.96 ± 0.03 


0.00001 ± 0.00005 


0.44 


0.28 


4/5 


100 


176-256 


-0.9 ±0.3 


0.0005 ± 0.0003 


0.67 


0.28 



TABLE II: The measured scaling exponents /i and v (with one-sigma statistical uncertainties) for 3-D numerical simulations of 
domain wall networks with different expansion rates (A) and wall thicknesses (Wo). The third column shows the range of the 
part of the simulation that was actually used in order to fit for the scaling exponent. The last two columns show the directly 
measured asymptotic values of ^c/?) and 711, which can be related to the macroscopic parameters of the analytic model. All 
simulations were done in 512"^ boxes with a — ?> and /3 = 0. Each value of the scaling exponent is obtained by averaging over 
25 simulations with different random initial conditions. 




FIG. 4: The evolution of the domain wall density (pu, = A/V) 
and velocity ((7^)^) as a function of conformal time, for the 
Wo = 100 boxes of Table HIl corresponding to different values 
of the expansion rate, A — 1/2,4/5. The plotted curves are 
the average of 25 different simulations with different (random) 
initial conditions. 



Bearing all this in mind, we believe that these re- 
sults are consistent with the notion that the networks 
should eventually reach the same attractor solution re- 
gardless of the wall thickness. An example can be seen 
in Fig. 131 which shows the results of simulations with 
Wq = 10, 50, 100 and all other parameters kept un- 
changed. 

In practical terms one is therefore justified in using the 
smallest wall thickness compatible with an accurate iden- 
tification of the walls. Previous results suggest |5j] that 
this is around Wq — 10, which is therefore the thickness 
used elsewhere in this paper. Nevertheless, similar re- 
sults could have been obtained with a larger thickness, 
albeit with some additional computational cost (compare 
Fig. Hwith Figs, mil and El . 



V. CALIBRATING THE VOS MODEL 

Table IIIII and Figs. [S] and [S] show the compared re- 
sults of very large 3D PRS simulations for three different 
cosmological expansion rates (parametrized by A): the 
usual radiation and matter eras plus and additional fast 
expansion era (with A = 4/5). As was the case with the 
simulations described in the previous section, we find no 
deviation from the scaling behavior. 

Moreover, comparing the scaling exponents for the cor- 
relation length in the 512^ and 1024^, and also for the 
a = 3, 256^ case of in Table U one can observe that the 
exponents tend to become closer to /i = — 1 as we increase 
the box size. This further supports our suggestion, in the 
previous section, that deviations from scaling are likely 
to be due to the limited range of the simulations. A fur- 
ther contributing factor may the the fact that the initial 
conditions that are normally used in these simulations 
produce networks that are quite far from scaling, and 
therefore the numerical evolution of these networks will 
require a considerable dynamic range before scaling can 
be reached. 

Given the large size and dynamic range of these simu- 
lations, in this table we now show the values of ^c/fl and 
"fv with error bars. As before we assume that the net- 
works are scaling exactly, and the values are calculated 
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Box Size 


A 


Fit range (77) 




u 




711 


512^ 


1/2 


38.5-256 


-0.99 ±0.05 


-0.0001 ± 0.0002 


0.60 ± 0.05 


0.46 ± 0.04 


1024^ 


1/2 


21-512 


-0.99 ± 0.05 


-0.0001 ± 0.0001 


0.64 ± 0.03 


0.48 ± 0.07 


512^ 


2/3 


26-256 


-0.97 ±0.04 


-0.0000 ± 0.0001 


0.54 ±0.01 


0.37 ±0.02 


1024^ 


2/3 


21-512 


-0.98 ±0.02 


-0.00001 ± 0.00005 


0.62 ±0.03 


0.37 ±0.03 


512^ 


4/5 


21-256 


-0.96 ±0.03 


0.00001 ± 0.00005 


0.44 ±0.01 


0.28 ±0.02 


1024^ 


4/5 


26-512 


-0.99 ±0.03 


0.00001 ± 0.00003 


0.50 ±0.02 


0.29 ±0.07 



TABLE III: The measured scaling properties of 3D PRS (a = 3, /? = 0) numerical simulations of domain wall networks with 
different box sizes and expansion rates (A). The third column shows the range of the part of the simulation that was actually 
used in order to fit for the scaling exponents /i and v. The last two columns show the directly measured asymptotic values of 
^c/j? and 7U, which can be related to the macroscopic parameters of the analytic model. All simulations have a constant wall 
thickness Wo = 10 and all error bars are one-sigma statistical uncertainties obtained by averaging over 25 simulations with 
different random initial conditions. 



X=\I2 - 




X=l/2 - 
>.=2/3 - 





X=2/3 ■ 



FIG. 5: The evolution of the domain wall density (p„ = ^/V) 
and velocity ((7«)^) as a function of conformal time, for the 
512^ boxes of Table HIl corresponding to different cosmological 
epochs: radiation era (A = 1/2), matter era (A = 2/3) and a 
faster expansion epoch with A = 4/5. The plotted curves are 
the average of 25 different simulations with different (random) 
initial conditions. 



FIG. 6: Same as Fig. [S] for the 1024^ simulations. 



from the last few timesteps in each simulation. We have 
experimented with several ways of calculating these val- 
ues, for example first averaging the 25 simulations and 
then calculating the scaling values from the averaged run 
versus calculating the scaling values from the individual 
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A 




kw 


1/2 


0.7 ±0.3 


0.8 ±0.1 


2/3 


0.5 ±0.2 


1.2 ±0.1 


4/5 


0.2 ±0.3 


1.6 ±0.2 


Weighted mean 


0.5 ±0.2 


1.1 ±0.3 



TABLE IV: The calculated values of the VOS model parame- 
ters and km for the 1024"^ boxes of Table Hill The final line 
show the best fit (weighted mean) result obtained by combin- 
ing the values from the three cosmological epochs. 



simulations and then determining the average of these 
values, and found that provided one has a large enough 
sample (typically involving the last 100 timesteps of each 
simulation) the results are remarkably consistent. The 
uncertainties shown are purely the one-sigma statistical 
errors for each set of 25 simulations. 

As expected, the results show that a faster expansion 
rate increases the damping on the walls, and therefore 
leads to smaller wall velocities and a larger wall density 
(which corresponds to a smaller correlation length). The 
agreement between the scaling values obtained from the 
512'^ and 1024'^ is also quite encouraging, and suggests 
that we are indeed seeing the network's attractor scaling 
solution. Nevertheless, it should be pointed out that the 
measurements of the velocities are in closer agreement 
than those of the correlation length. This is likely to 
be due to the fact that our algorithm for identifying the 
walls and adding their areas is not completely accurate. 

These results can now be used to calibrate the VOS 
walls model. Recall that in this model the scaling solu- 
tion is parametrized by the expansion rate A (such that 
a oc t^) and the phenomenological parameters and . 
Given these parameters, the predicted values for the scal- 
ing arameters e and v are given by Eqs. pTHT^ . Using 
the results of Table Hn] we can trivially obtain the value 
of V, and e is also easily obtained from 



(1 - X)v ' 



(15) 



from these one finally obtains the numerically measured 
values of and km. The results of this analysis, using 
the data from the 1024^ runs, are shown in Table IIVI 
We first calculated these values separately for each of 
the three cosmological epochs that we simulated, and we 
finally combined the three into a final, weighted mean set 
of results. Since these purely statistical errors may be 
underestimates, we rescaled the variance in the standard 
way, by multiplying it by the chi-squared per degree of 
freedom. Our final calibrated model parameters are 



and 



= 0.5 ±0.2 



1.1 ±0.3; 



(16) 



(17) 



these results are remarkably consistent with the previous, 
more qualitative analysis in [18], which had found ~ 
0.5 and k^ ^ 0.9. 



VI. CONCLUSIONS 

We took advantage of recent improvements in comput- 
ing power to numerically study the evolution of the sim- 
plest domain wall networks. Having carried out sets high- 
resolution, large dynamic range simulations with various 
amounts of damping, we find strong support for the sug- 
gestion that the attractor solution for the evolution of 
these networks is a linear scaling solution, with ^phys oc t 
(or equivalently v) ^-nd v = const. Our results sug- 
gest that previous hints of deviations from this behavior 
may have been due to the limited dynamical range of 
those simulations. 

Moreover, we have used the results of the largest 
(1024^^) of our simulations to provide a calibration for the 
velocity-dependent one-scale model for domain walls. As 
a consistency check, we have used results from simula- 
tions in three different cosmological epochs, even though 
the model only has two free parameters, one of which 
quantifies the network's energy loss rate while the other 
describes the (curvature-related) forces acting on the 
walls. Given the conceptual simplicity of the analytic 
model, we believe that the present numerical results sup- 
port its validity, and suggest that it can be reliably used 
as a tool to study the cosmological consequences of these 
networks in quantitative detail. 

This combination of analytical and numerical tech- 
niques, leading to a detailed calibration of a model 
(which so far had only been carried out for cosmic strings 
|2ll [26j ) can in principle be extended to networks of do- 
main walls with junctions. The larger number of degrees 
of freedom (corresponding to additional scalar fields) 
makes the study of these models numerically trickier 
since in most cases the factor limiting the size of the boxes 
that can be simulated is memory rather than time, but 
otherwise our methods are directly appicable there. An- 
other case of interest is that of semilocal string networks: 
here some steps towards an accurate model calibration 
have been taken recently [13], and we shall return to it 
in a subsequent publication. 
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